VAST MULTIPLICITY OF VERY SINGULAR SELF-SIMILAR 
SOLUTIONS OF A SEMILINEAR HIGHER-ORDER DIFFUSION 
EQUATION WITH TIME-DEPENDENT ABSORPTION 



V.A. GALAKTIONOV 

Abstract. As a basic model, the Cauchy problem in X M , for the 2mth-order 
semilinear parabolic equation of the diffusion-absorption type 

ut = -(-A)"u - eiulP-^u, with p > 1, a > 0, m > 2, 

with singular initial data uo ^ such that uq{x) — for any x 7^ is studied. The 
additional multiplier h{t) — t" as t ^ in the absorption term plays a role of time- 
dependent non-homogeneous potential that affects the strength of the absorption term 
in the PDE. Existence and nonexistence of the corresponding very singular solutions 
(VSSs) is studied. For m = 1 and h{t) = 1, first nonexistence result for p > po = 1 + 
was proved in the celebrated paper by Brezis and Friedman in 1983. Existence of VSSs 
in the complement interval 1 < p < pa was established in the middle of the 1980s. 

The main goal is to justify that, in the subcritical range l<p<po = l+ ^ 
there exists a finite number of different VSSs of the self-similar form 

u,{x,t)=t-PV{y), y = x/t^, P = 

where each V is an exponentially decaying as y cxd solution of the elliptic equation 

_(_A)™-l/+_l.y.VX^-|-/3y-|l/|P-V = in M^. 

Complicated families of VSSs in ID and also non-radial VSS patterns in are detected. 
Some of these VSS profiles VJ are shown to bifurcate from at the bifurcation exponents 

P; = l + 2^^^, where Z = 0,1, 2,.... 



1. Introduction: VSSs for higher-order diffusion-absorption PDEs 

1.1. Diffusion-absorption model. Our basic model is the 2mth-order semilinear heat 
equation in x with non-autonomous (non-homogeneous) absorption term 

(1.1) ut = -{-/^Tu-h{t)\u\^~\ (p>l, m>2), 

where the given function h{t) satisfies 

hit) > 0, h{t) ^0 as t ^ 0. 



Date: January 27, 2009. 

1991 Mathematics Subject Classification. 35K55, 35K40, 35K65. 

Key words and phrases. The Cauchy problem, diffusion equations with absorption, initial Dirac mass, 
very singular solutions, existence, nonexistence, bifurcations, branching. 

1 



For m = 1 and h = 1, the critical Fujita-like exponent po = 1 + jj and existence- 
nonexistence of very singular solutions (VSSs) was derived since the 1980s; see first prin- 
cipal results of the 1980s in H [151 El [121 [2ll [221 [23], and a long list of references in 
[T8l Ch. 4]. In the higher-order case m > 2 and h = 1, such an investigation has been 
performed in [IH [IZl [IS] . The influence to VSSs of non-homogeneous potentials h(t) was 
first studied in Marcus- Veron [26] in the case m = 1; see more recent extensions in |31j . 



1.2. Similarity VSSs: towards a non-variational problem. The layout of the paper 
is as follows. We investigate the phenomena associated with the non-homogeneous term 
and consider the power case 

(1.2) h(t) = t", with the fixed exponent a > —1 (so that j^hif) dt < oo). 

In Section [21 for the higher-order equation (11. ip with m > 2, we study the existence and 
multiplicity of similarity solutions and show that, in the subcritical range p G (l,po)) with 
the critical (Fujita-like) exponent 

_ . 2m(l+«) 

Po — -L H ]v — ' 

there exist VSSs of the similarity form 

(1.3) u,{x,t)=t-^V{y), y = x/t^, with /? = i±f > 0, 
where is a non-trivial radial solution of the elliptic equation 



Bil/ - \V\'P-^V = -(-A)'"1/ + ^yW + 13V - \V\-p~^V = in 



1.4) 

V{y) decays exponentially fast as \y\ ^ oo. 
Passing to the limit t ^ 0^ in (11.31) yields the initial function satisfying 

(1.5) u(x,0) = for all x 0, 

which nevertheless is not a finite measure, and this reinforces the notion of very singular 
solutions for (II. 3p . Condition in (ll.4p on exponential decay at infinity is naturally enforced 
by introducing weighted and Sobolev spaces; details are given in in Section [H 

1.3. On main results. The main goal is to detect the whole range of VSS profiles by 
using spectral properties of the linear operator Bi in (II. 4p and bifurcation theory. This 
determines a countable number of p-bifurcation branches of VSS profiles in the subcritical 
range p < Po that, for any a > — 1, are originated at the bifurcation points 

(1.6) pi = l+'-^, / = 0,1,2,.... 
This also yields a-bifurcation points, 

(1.7) = (£z]M±0 _ 1, / = o,l,2,.... 

In Section [3l using numerical experiments, we show that, surprisingly, the whole picture 
of branches of VSS similarity profiles is essentially more complicated than that for a = 
studied in [19J. It turns out that even in ID for sufficiently large a > 0, the whole set 
of VSS profiles consists of various families of solutions of different geometric shapes. In 
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addition, this includes new branches that are expected to appear at saddle-node (turning) 
bifurcation points in a. In Section HJ for future use for non-power potentials such as 



(1.8) /i,(t)=e 



we are particularly interested in the study of the limit behaviour of such VSSs as a +00. 
Some nonexistence results of VSSs in the range p > po are proved in Section [51 

It is necessary to mention that for, any m > 2, the problem (11.41) even in ID is not 
variational (it is for m = 1 only), so the operators there are not potential in any topology. 
Moreover, even the linear principal part Bi therein in not self-adjoint (symmetric) in any 
weighted space L^. In addition, the corresponding parabolic flow is not order-preserving 
in view of the lack of the Maximum Principle (again applies for m = 1 only). 

Therefore, for describing the complicated discrete sets {K^j.}, where cxfc is a multiindex 
(see Section |3]), we are going to use a machinery of various bifurcation- branching methods 
of nonlinear analysis. Nevertheless, many of our final conclusions on the global behaviour 
of bifurcation diagrams of V^'s remain formal, and we do not expect that a reasonably 
simple justification can be achieved soon. Therefore, we heavily rely on numerical con- 
struction of, at least, am couple of hundreds of similarity profiles. These numerics are 
also not that easy at all. Often, the reliable numerics will be the only tool to detect this 
vast multiplicity of VSS profiles of (II. 4p far away from bifurcation points. 

2. VSSs PROFILES: LOCAL BIFURCATION THEORY 
Thus, consider the Cauchy problems for the PDE (II. ip with the power potential (II. ip . 

2.1. The fundamental solution of the poly-harmonic equation. This is 

N 1 

(2.1) b{x,t)=t-^F{y), y = x/t^, 
being the fundamental solution of the linear poly-harmonic equation 

(2.2) ut = -{-Aru. 

The rescaled kernel F is then the unique radial solution of the elliptic equation 

(2.3) BF = -i-ArF + CoF = inm^, f F = 1; C, = ^y.V + £l, 
satisfying for some positive constants D > 1 and d > depending on m and p!0] 

(2.4) \F{y)\<DF{y)^Du,e-'\y\'' in R"", a = e (1,2), 

Ui > being a normalization constant such that j F = 1. 

It turns out that the linear operator Bi in equation (II. 4p is connected with the operator 
(12. 3p for the rescaled kernel F as follows: 

(2.5) Bi = B + c^I, where c, = 
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2.2. The point spectrum of the non self-adjoint operator B. We consider the 
hnear operator B given in fl2.3p in the weighted space L^(R^) with the exponentiaUy 
growing weight function 

(2.6) p{y) = e'^l^l" > in M^, 

where a G {0,2d) is a sufficiently small constant. We ascribe to B the domain Hp"^(R,^) 
being a Hilbert space with the norm 



\v 



''-Ip{y)ElZo\D'viy)\'dy, 



induced by the corresponding inner product. We have Hp"^ C C L^. The spectral 
properties B follows [9]: 

Lemma 2.1. (i) B : H^"^ —>■ L"^ is a bounded linear operator with the real point spectrum 

(2.7) a(B) = {A, = -i, / = 0,1,2,...}. 
The eigenvalues A/ have finite multiplicity with eigenf unctions 

(2.8) My) = ^ D^F{y), with any = / > 0. 

(ii) The set of eigenf unctions $ = {ipp, \[3\ = 0, 1, 2, ...} is complete and closed in L^. 
In the classical second-order case m = 1, we have that 

F{y) = (47r)~fe-3lJ^I' 
is the rescaled positive Gaussian kernel and the eigenfunctions are 

My)=^~''^'^'Hf,iy), 

where Hj3 are separable Hermite polynomials in [3l p. 48]. The operator B, with 



the domain H"^ and the weig ht p = el^^l /^ is self-adjoint and the eigenfunctions form an 



orthogonal basis in L^. 



Lemma [271] gives the centre and stable subspaces of B, E'^ = Span{-?/'o = F}, = 
SpanfT/-^, \(3\ > 0}. 

2.3. The polynomial eigenfunctions of the adjoint operator B*. Consider the 
adjoint operator to B, 

(2.9) B* = -(-A)--^yV. 

For m = 1, the following representation holds: 

B* = ^ V ■ (p*V), with V{B*) = H^,, where p*{y) = e"! 1^1", 

is self-adjoint in L^* and has a discrete spectrum. The eigenfunctions form an orthonormal 
basis in L^* and the classical Hilbert-Schmidt theory applies |3j. 

For m > 1, we consider B* in L^, with the exponentially decaying weight function 

P*iy) = ±^^e-^\y\'' >o. 
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Lemma 2.2. (i) B* : H'^P- —>■ L^, is a bounded linear operator with the same spectrum as 
B, ( |i^.7p . The eigenf unctions ip*p{y) with \(3\ = I are Ith-order polynomials 

(2.10) rp{y) = ^.V + Y}M^ jA-^T'y^] ■ 

(ii) The set {ip*^} is complete in L^*. 

From this definition of the adjoint eigenfunctions, the orthonormahty condition holds 

(2.11) {^p,r,) = 5p^, 

where (■, ■) denotes the standard inner product. For m = 1, both (12. 7p and (I2.10p are 
well-known properties of the separable Hermite polynomials, [3]. 

2.4. Bifurcations at p = pf. local existence of the VSSs. It follows from (12. 5p and 
the above lemmas that the only possible bifurcation points of VSSs are obtained from 

(2.12) ci = -A, =^ p = p^ = i + Mlj!), / = o,l,2,.... 

We justify these bifurcation phenomena. Taking p near the critical values as defined in 
(12.121) . we look for small solutions of the problem (II. 4p . At p = pi, the linear operator Bi 
has a nontrivial kernel, hence, the following result: 

Proposition 2.1. Let for an integer I > 0, the eigenvalue A/ = — 2~; of operator { \2.!^ he 
of odd multiplicity. Then the exponent in {\2.1^ is a bifurcation point for ( |i.^p . 

Proof. This result is standard in bifurcation theory; see [SJ p. 381] and similar results in 
a related problem in [191 § 6] . Necessary spectral properties of the linearized operator are 
given in [H]. We present some comments concerning the linearized operators involved. 
Consider in the equation 

(2.13) BV = -{l+Ci)V +\V\P-^V, B = Bi-(l+ci)/ = B-/. 

The spectrum of B is a translation of that of B, cr(B) = { — 1 ^ 2m }' ^"^^ consists of strictly 
negative eigenvalues. The inverse operator B~^ is known to be compact (Proposition 2.4 
in [9]). Therefore, in the corresponding integral equation 

(2.14) V = A{V) = -(1 + ci)B- V + B-^ll^l^'- V, 

the right-hand side contains a compact Hammerstein operator in L9(M^) space for some 
g > 1 [211 p. 38] (see details on the resolvent of B in [6]). In this application, there exist 
certain technical difficulties in checking compactness of this Hammerstein operators in 
weighted L'^-spaces over whole M^. As an alternative, we can use Ladyzhenskii's theorem 
[211 p. 34] establishing compactness in C. 

To avoid this technicalities, we use [8, Thm. 28.1] where no assumptions on compactness 
of the vector field are necessary. The main hypothesis therein is oriented to the linearized 
operator in (12.140 that is assumed to be Fredholm of index zero at bifurcation values 
(which is true by [9|). Then the result is first obtained for truncated uniformly bounded 
nonlinearity instead of |V^|''~^V^ (see [l9l p. 1088]), for which all the hypotheses are valid. 
Then bifurcations in the truncated problem (12.140 are always guaranteed if the derivative 



A'(0) = —(1 + ci)B~^ has the eigenvalue 1 of odd multiphcity. (cf. also ^24, p. 196] for 
compact integral operators). 

Thus, bifurcations in the problem (12.141) occur if the derivative A'(0) = —(1 + ci)B~^ 
has the eigenvalue 1 of odd multiphcity. Since cr(A'(0)) = {(1 + ci)/(l + 2^)}, we arrive 
at critical values fl2.12p . By construction, the solutions of (12.141) for p ^ pi are small in 
and, as can be seen from the properties of the inverse operator, in if^'". Since the 
weight p{y) is a monotone growing function as \y\ 00, this implies that V G H^"^ in 
the subcritical Sobolev range 

(2-15) l<p<ps = ^ 

(which is true in our VSS range p < po) is a uniformly bounded, continuous function by 
standard elliptic regularity results and embedding theorems; these can be found in [27] 
and ^32j. For < 2m, the result is straightforward in view of embedding [32l p. 5] 

i7™(M^) C C(M^) {N < 2m). 

Let us mention related boundedness results of parabolic orbits of (11.41) {a = 0) in the 
range (12.151) obtained in fi9[ § 2] by Henry's version of Gronwall's inequality with power 
kernel and scaling arguments. 

Therefore, for p ^ pi, we have bounded, uniformly small solutions of (12.141) only. By 
interior regularity results for elliptic equations, these solutions are smooth enough to be 
classical ones of the differential equation (12.131) . □ 

Thus, / = is always a bifurcation point since po = is simple. In general, for 
I = 1,2,..., the odd multiplicity occurs depending on the dimension A^. In particular, for 
1 = 1, the multiplicity is N, and for / = 2, it is i^M+ll, the case of the even multiplicity 
of A/, an extra analysis is necessary to guarantee that a bifurcation occurs, [25j. It is 
important that, for key applications, namely, for = 1 and for the radial setting in M^, 
the eigenvalues (12.71) are simple and (12.121) are bifurcation points. 

Since the nonlinear perturbation term in the integral equation (12.141) is an odd suf- 
ficiently smooth operator, we easily obtain the following result describing the local be- 
haviour of bifurcation branches, see EH and EHl Ch. 81. 



Proposition 2.2. Let Xi be a simple eigenvalue o/B with eigenfunction ipi. Denoting 

(2.16) Ki = mp-%,^t), 



we have that problem ( |i.-^p has (i) precisely two small solutions forp ^ pj and no solutions 
for p p^ if Ki > 0, and (ii) precisely two small solutions for p ^ p'l and no solutions 
for p ^ p'[ if Ki < 0. 

In order to describe the asymptotics of solutions as p — >• pz, we apply the Lyapunov- 
Schmidt method ([251 Ch. 8]) to equation (12.141) with the operator A being different iable at 
0. Since under the assumptions of Proposition 12.21 the kernel Eq = ker A'(0) = Span {ipi} 
is one-dimensional, denoting by Ei the complementary (orthogonal to tp'l) invariant sub- 
space, we set V = Vq + Vi, where Vq = eiipi G Eq and 

^1 = Y^k^i^ki^k e El. 



Let Pq and Pi, Pq + Pi = I, be projections onto Eq and Ei respectively. Projecting (12.141) 
with n = oo onto Eq yields 

(2.17) ^,e, = {B-\\Vr^VUt), 1i = 1-^ = -jtt£^^, 

where s = pi — p. By the general bifurcation theory (see e.g. [23 p. 355], [8, p. 383], and 
branching approaches in [33]; note that operator A'(0) is Fredholm of index zero), the 
equation for Vi can be solved and this gives Vi = o{ei) as ei — > 0, so that ei is calculated 
from the Lyapunov bifurcation equation (12.171) as follows 

7/6/ = \eirhi{B-'\^Pir%,^!) + o{\ei\n =^ \eir' = ci[ipi-p)+o{l)], 



where q = ^Ji^jt^l^-,^^ . Here, we have performed calculations as follows: 

(B-vf , = (^f , (B*)- vn = + h)- 

It is natural to require ni > 0, though this is not always true for a > 0; see a closed curve 
of solutions in Figure [2] (for a = the branches were detected to be always monotone, 
[TU]). In view of the orthonormality property (12. lip , for p = 1 we have = 1, so that by 
continuity we can guarantee that 

(2.18) Hi > at least for all p ^ 1+. 

Thus, we obtain a countable sequence of bifurcation points (12.121) satisfying — > 1"*" as / — > 
oo, with typical pitch- fork bifurcation branches appearing in a left-hand neighbourhood, 
for p < pi. The behaviour of solutions in H"^"^ and uniformly takes the form 

(2.19) Vi{y) = ± [q{pi -p)f^^-'^ (My) + o(l)) as p ^ p^. 

We now prove the main result concerning "local" existence and stability of the VSS 
solution with the similarity profile Vo{y) corresponding to the first bifurcation point, 
p = Pq. If kq > 0, as expected, then two bifurcation branches exist for p < pq. 

Theorem 2.1. For p ^ p^, the problem ( [i.^/p admits a solution Vo(?/) ^ provided that 
^ is small enough, and then it is an asymptotically stable stationary solution. 

Proof. As we have shown, a continuous branch bifurcating at p = exists if 

(2.20) Ko = (|^or-Vo,V'S) = j\F\'^F > (^0* = !)• 

In view of the positivity dominance of the rescaled fundamental solution F, J F = 1, we 
have that (I2.20p holds by continuity provided that ^ -C 1. Therefore, in this case there 
exists a solution (12.19^ with / = satisfying for small s = po — p > uniformly 

(2.21) Voiy) = icos)^^[Fiy)+oil)], co = 

We now estimate the spectrum of the corresponding linearized operator 

(2.22) Do = Bi-p|\/or'/. 

Some of the eigenvalues of the operator (12.220 follow from the original PDF (II. ip . For 
instance, the stable eigenspace with \ = —1, ip = -^^ziVq + ^y ■ VVq G L^, follows from 



the time-translational invariance of the PDE. For = 1, translations in x yield another 
pair A = —3^, = V^y G L^. For > 1, in the non-radial setting, this p has multiplicity 

with eigenf unctions Voy-. These are not the first pair with the maximal real part. 

Bearing in mind that the spectrum of the unperturbed operator B is real, (12. 7p . and 
has the unique, non-hyperbolic eigenvalue Aq = 0, we use (I2.2ip to obtain 

(2.23) Do = B + s(l + o(l))C, 

where, as it follows from (I2.20p and (12.211) aX p = po, the perturbation has the form 
(2-24) C^^^(1-^^\F\^)I. 
Therefore, we consider the spectrum of the perturbed operator 

(2.25) Do = B + sC. 
Since (B — /)^^C is bounded, 

(Do - I)-' = [I + s{B - I)-'C]-\B - I)-' 

is compact for small |s| as the product of a compact and bounded operators. Hence, 
Do also has only a discrete spectrum. By the classical perturbation theory of linear 
operators (see e.g. [20]), the eigenvalues and eigenvectors of Dq can be constructed as a 
perturbation of the discrete spectrum (t(B) consisting of eigenvalues of finite multiplicity. 
We are interested in the perturbation of the first simple eigenvalue po = 0. Setting 

^0 = sfJ'0 + o{s), -ipo = ipo + sifo + o{s) as s 

and substituting these expansions in the eigenvalue equation Uoipo = ^oi^o yields 

(2.26) B(^o = (-C + /io/)V'o- 
We then obtain the solvability (orthogonality) condition 

{{-C + i^oI)^o,ro) = =^ /xo = (CF,l). 

Using (12:241) yields 

f^o = -£< 0. 

Therefore, RcAq < — ^ < for all p ^ Pq ■ Since, with these properties of the spectrum, 
the perturbation (I2.22p of B remains a sectorial operator with o"(Do) C {Re A < — ^} 
and ||e^i^||£ < Ce-^^P^-P^^/^"" in the norm of C{Hl'^,Hl'^) [13], V^iy) is exponentially 
stable in Hl"^. □ 

We expect that condition ( 12.201) remains valid for any m and so that Vo(?/) is stable 
without the restriction 2m ^ A^. We have a numerical support for this, but, as yet, no rig- 
orous proof exists. Possibly, to check conditions such as (12.201) we must currently rely on 
numerical evidence and then, as often happens in spectral theory and applications. The- 
orem 12.11 can be established with a hybrid analytic-computational proof. We also expect 
that the whole branch bifurcating from p = po remains stable for all p G {l,po), though 
the proof would require to establish that the discrete spectrum o"(Do) never touches the 
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imaginary axis. In particular, this open problem means that a new (nonlinear) saddle- 
node bifurcation never occurs on this po-branch, i.e., it does not have turning points. 
For the variational problem with m = 1, this is valid [30] as well as for ordinary differ- 
ential higher-order equations with self-adjoint positive operators of special structure of 
quasi-derivatives [HES]; see also properties of bifurcation branches in Berger [21 p. 380]. 

Further, it is easy to see that the other bifurcation branches are unstable. Taking any 
/ > 1, instead of (12.231) we now have 

= Bi - p\Vir'l = B+[c,- spiciil^ir' + 0(1))] J, s = pi-p. 

From the definition of Bi, (12. Sp . Ci > for all p ^ pi, thus Vi for / > 1 is unstable. 

2.5. Remark on global bifurcation diagrams. Global bifurcation results concerning 
continuous branches of solutions were already given in Krasnosel'skii [2ll p. 196] (the 
first Russian edition was published in 1956). Concerning further results and extensions, 
see references in [HI Ch. 10] (especially, see [H p. 401] for typical global continuation of 
bifurcation branches), and also [25^, § 56.4]. These approaches deal with integral equations 
with compact operators such as (12.141) . 

In the present non- variational problem, the main open problem of concern is to establish 
under which conditions the p-branch originated at some p = pi can be extended up to 
p = 1~, so that cannot end up at another bifurcation point pk < pi- Actually, this can 
happen for a > 0; see Figure [21 where the closed branch is originated aX p = po and ends 
up at p = p2- 



(3.1 



3. p- AND a-BIFURCATIONS AND BRANCHES: VSS MULTIPLICITY 

We consider the fourth-order equation (II. 4p in ID with m = 2, i.e., the ODE problem 
-V^^'^ + lyV' + f3V -\V\P-^V = in R (/?=i±f), 
V{y) has exponential decay as ?/ oo. 



Note that the condition of the exponential decay is crucial for VSS setting, since the ODE 
in (13. ip admits continuous families of solutions with algebraic power decay such as 



2m 



(3.2) V{y) = C±\yr—{l + o{l)) as y ^ ±oo, 

where C± G M are, in general, arbitrary constants; see an example in [T9l § 8] for the case 
a = 0. In numerical applications, the exponentially decaying solutions of (13. ip are always 
clearly oscillatory for ?/ ^ 1 that strongly differ them from those satisfying (13. 2p . 

For even profiles Vq, V2,... satisfying (13. ip that we pay more attention to, we impose 



the symmetry conditions at the origin 

(3.3) V^'(O) = V"'{0) = 0. 

For the odd profiles Vi, V3,... , we pose the anti-symmetry conditions 

(3.4) V{0) = V"{0) = 0. 



p^-branch for a^O, m=2 Deformation of V^(y) profiles for a=0, m=2 




p y 

(a) po-branch (b) Vq profiles 

Figure 1 . The first po-branch of solutions of (13. ip for m = 2 and a = 0; the 
branch (a) and deformation of profiles Vo{y) with p, (b). 



3.1. p-branches for fixed a > 0. As we have mentioned, the autonomous case a = 
is rather well-understood, [19] . For example, in Figure [1], we present the first strictly 
monotone po-branch of similarity profiles Vo{y) for a = 0. This branch blows up as 
p — > 1~ according to the asymptotics in p[9| p. 1091]. 

For a > 0, the p-bifurcation branches can essentially change their topology and can be 
closed curves, as Figure [2] demonstrates for a = 1. Here, the first pg-branch is initiated 
at po = 9 and has the saddle-node (turning) point at 

Ps-n = 3.51... . 

Eventually, the branch ends up at the third critical point (11.61) . i.e., at 

P2 = 1 + I = 3.666.... 

For convenience, in Figure [3l we give the enlarged structure of the almost vertical part of 
the branch in Figure [2]^a) and VSS profiles on it for p G [Ps-n,P2]- 

3.2. On a-branches for fixed p. A typical example is presented in Figure HI where we 
show the first a-branch for N = 1, p = 2, that, according to (II. 7p . is originated at the 
bifurcation point 

«o = -i (/ = 0). 

We observe from (b) that the profiles get wider as a increases. This is a natural phenom- 
enon that will play an important role in what follows. It should be noted that Figure (b) 
does not explain jump-discontinuities that are not only possible but often happen in such 
numerics (in view of huge density of various branches of solutions), and were observed 
even with the enhanced Tolerances 10~^ and the step Aa = 10~^. In fact, (a) shows for 
a 1 a certain non-smoothness of the branch that we cannot explain. To be honest, 
we cannot guarantee that here this is not associated with a discontinuous jump to other 
neghbouring branches. The a-branching deserves both extra numerical and analytical 
study. 
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Closed Pj^-branch, m-2, a-^ 



Profiles o1 P|^-branch from p^^ to p^_^, m-2, a=1 




p .1+8/3=3.6 




10 12 



(a) po-branch (b) Vq profiles 

Figure 2. The closed po-branch of solutions of ()3.ip for m = 2, and a = 1; the 
branch (a) and deformation of profiles Vo{y), (b). 



p,-brancfi from s-n bifurcation to p m-2, a-l 



Profiles of p -branch for p e (3.51 ,p ), m-2, a-l 




(a) po-branch 
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(b) Vq profiles 



Figure 3. The po-branch of solutions of (j3.ip between two bifurcations; p G 
[3.51,3.66], m = 2, a = 1; the branch (a) and deformation of profiles Vo(y), (b). 



3.3. Further bifurcations and discussion. We denote the basic set of VSS profiles by 
(3.5) basic spectrum: Vq, Vi, V2, V^, ... . 

Each Vk{y) satisfies the approximate Sturmian property, i.e., has precisely k dominant 
extrema. For m = 1, this is true sharply, by the Maximum Principle. For m > 2, by 
"dominant" we mean those minima or maxima that are different from smaller ones in 
the oscillatory tail of the solutions. This approximate Sturm's zero (extremum) property 
for iV = 1 is associated with the fact that in the ODE (11.41) with m > 2, the leading 
2mth-order operator is a multiplication of m positive second-order ones, 
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for which Sturm's property is true concerning its eigenf unctions; see general conclusions 
in EUias pjj, also applications to bifurcation diagrams of interest here in P,|29]. It turns 
out that lower-order perturbations in (11.41) create only sufficiently small extrema and 
infinitely many zeros in the exponential tails of solutions. 

For a > 0, the basic profiles are essentially deformed. For instance, a few of such 
profiles (already available in Figure |l](b)) are shown in Figure [5] for various a > 0. 

A principal new feature of the problem with the exponent a is that there may occur 
other saddle-node bifurcations in the ODE in (13. ip . The origin of extra a-bifurcations 
can be seen as follows. The leading operator in (11.41) . 

is coercive and negative in the topology of the Sobolev space H'^"^{Bi) fl L^^^{Bi) in the 
subcritical range (12.151) (this is enough for = 1). Here Bi is the ball of the radius I > 
in R^. It follows from Lemma [2.1( i) that a similar result remains true in the weighted 
space L^(]R^), and moreover, the first linear term ■ W can be estimated via the 
leading operator. We then observe that the crucial part is played by the second term 

The multiplier gets arbitrarily large positive for a ^ 1, and hence can produce a 
bifurcation, where the negative and positive operators are sufficiently balanced; see more 
precise estimates below. 

Let us discuss THE numerical evidence more systematically. In Figure O we show the 
first VSS profiles Vo{y) satisfying (ll.4p for p = 3, = 1, and for a number of different 
values of a. The bold line corresponds to Vq in the classic case a = studied in [19]. 

In Figure [71 we show twelve different even VSS profiles V2k{y) for p = 2, N = 1, and 
a = 4. Note that according to critical exponents (11.71) . there exist precisely ten critical 
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V/ a=7 
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y 

Figure 5. Various VSS profiles associated with the first basic pattern Vo(y) 
for p = 2, = 1, and different values of a € [0, 7]. 

Deformation of profiles Vjj(y) with a, m=2, p=3 




y 

Figure 6. Deformation with a of the first profile Vo{y) of (j3.ip . with p = 3, 
for a G [-0.5,7]. 

values below 4, at which such profiles can be originated, 

"O = — tt2 = — ^4 = \, ■■■ , ai6 = X' '^18 ^ T' 

On the other hand, according to (11.61) . there exist the same number of critical exponents 
pi that are above 2, 

po = 21, p2 = f , ... , Pie = pi8 = ff. 

It follows that at least two profiles in Figure [7] are not generated by standard p- or a-curves 
originated at bifurcation points. This mystery remains an open problem and needs extra 
analysis. We expect that some pairs of these profiles could be originated at saddle-node 
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Profiles for m=2, p=2, a=4 
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Figure 7. Twelve even VSS profiles V2k{y) for p = 2, iV = 1, and q = 4. 

Basic spectrum, p=1 .5, N=1 , a=0 
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V(y)^ 
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-1.5 
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Figure 8. Three VSS profiles Vq, Vi, and V2 for a = 0;p = 1.5 and N = 1. 

bifurcations at some cts-n > 0. It is seen that the profiles slightly "oscillate" about the 
constant equilibrium of the ODE in (13.11) which is 

V+ = (i±f ) ^ = 5 for p = 2 and a = 4. 

For convenience, below we present the standard multiplicity result for the autonomous 
case. In Figure [H we show first three solutions Vo{y), Vi{y), and V2{y) of the problem 
(13. ip with p = 1.5 and a = 0. As was shown in [19], these p-bifurcation branches exhaust 
the whole set of VSS profiles in the case a = 0. Notice that the mutual geometry of VSS 
profiles in the case a = essentially differs from that in Figure [71 
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Basic profiles V (y), m=2, p=1 .5, N=1 , a = 7 
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Figure 9. Several VSS profiles V2k{y) from satisfying (13. ip for p = 1.5 and a = 7. 



"Dipole" profiles ^/^^^^(y), m=2, p=1.5, N=1, a=7 
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-100' • • • • • • ' 
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Figure 10. Several odd VSS profiles V2k+i{y) satisfying (f3TT]) for p = 1.5 and 
a = 7. 

Let us more clearly show that, for a > 0, there occur other bifurcation-branching 
phenomena. In Figures [9] and [TOl we again consider the case p = 1.5, N = 1, and a = 7. 
For this sufficiently large a, according to ( 11. 7p . we observe several profiles from the basic 
family {V2fc} of even functions that are now concentrated about the constant equilibrium. 
Namely, in Figure [HI we show nine profiles {V2A:}. Figure [TUl demonstrates thirteen odd 
profiles from the family {V2fc+i}. In this case, instead of symmetry conditions (13.31) . we 
take the anti-symmetry ones (13. 4p . Then the solutions V{y) are continued in the odd 
manner, V{—y) = ~V{—y) for y < 0. 

Before performing some easy preliminary estimates, we show in Figure [TT] the case p = 4 
and a = 10, where a similar set of seven VSS patterns from the even family {V2fc} are 
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Two families of profiles, p=4, N=1 ,a=10 




y 

Figure 11. VSS profiles from {V2k} and {F^J for p = 4, iV = 1, and a = 10. 



Profiles V (y), m=2, p=2, a=10 




y 

Figure 12. Profiles from the family {V2fc} for m = p = 2, N = 1, and a = 10. 
presented. In addition, we see here another new family 

(3.6) {^o-i}) with multiindex ai = {+0, 1, — A;}, 

where 0, 1, and k in ai stand for an "effective number of intersections" of the profiles with 
three consecutive equilibria V+,0, and — V+. We will discuss this new family later on. 

In Figure fT2l we show the basic profiles V2k{y) and others for m = p = 2 and a = 10. 
This computation is important to treat the results in Figure H] on the a-branches. 

For future convenience, some of the results are natural to present for the rescaled 
equation by using the scaling 

(3.7) V^CV, y^cy, where 0^'^ = f3, c = C'^ = (3^ , 
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Profiles Vjj, m=2, p=1 .5, a.7 Profiles , m.2, p=1 .5, a=7 




10 20 30 40 50 60 70 10 20 30 40 50 60 



y y 

(a) profiles V2k (b) profiles 

Figure 13. Two different families of VSS profiles for m = 2, p = 1.5, = 1, 
and a = 7. 

that gives the following rescaled ODE: 

(3.8) (_i)™+V(2m) ^ iA_yV' + V ~ \V\P-^V = 0. 

The constant equilibria are now fixed and are independent of the parameters, 

(3.9) V± = ±1 (and = 0), 

so we do not observe such huge VSS profiles as, for instance, in Figured The correspond- 
ing rescaled families {V2k} and {K-j} for m = 2, p = 1.5, and a = 7 are shown in Figure 
[T3l(a) and (b). Figure shows another more complicated family of even VSS profile 

V^M, where ^2 = {+1, 1, -2, 1, +A;}. 

By boldface dotted line we denote therein a different profiles from the family {K-g}, with 
a distinct multiindex = {— 0, 1, +2, 1, — A;}. It seems that there are other families 
of VSS profiles with more involved geometric structures that are very difficult to catch 
numerically, so we stop the discussion at this moment. 

We expect that all the new profiles from families {K-i}, {K-a}' ^tc, are originated at 
saddle-node bifurcations at some a^^ > and cannot be extended to a = 0. Recall that 
such VSS solutions were not detected in the autonomous case a = 0, [19]. Numerics 
with Aa = 10~^ show some evidence concerning this though are also rather unstable 
demonstrating typical jumps between many neighbouring a-branches. We believe that 
such stable jumps definitely show existence of saddle- node points that are difficult to 
detect numerically. 

On countable sets via shooting arguments. As an illustration, consider the ODE 
in (13.11) (i.e., for m = 2 and = 1) in the symmetric case (13. 3p . so we study V2k and 
other even profiles. The shooting argument justifies existence of a connection of the 2D 
bundle near the origin, which is characterized by two symmetry conditions (13.31) . so that 
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VSS family of profiles V^^ ^ ^ ^ ^^{y), m=2, p=1 .5, a=40 
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Figure 14. Some profiles from the family {V^al m = 2, p = 1.5, and a = 40. 
the bundle has two parameters 

(3.10) /(O) and /"(O). 

We then need to intersect this bundle with the exponential bundle at y = +oo, where 
the linearized ODE admits the following WKBJ-type asymptotics (we omit power-like 
multipliers that are not important at this stage): 

V^'^ = + - =^ Viy) ~ e'^''''' with = 

Therefore, there exists a ID unstable bundle with positive = | 4~^/^, and a 2D stable 
oscillatory bundle with 

a± = ao(-i±i^) (Rea±<0). 

Thus, we observe here a typical difficult two-parameter shooting problem of 2D— s>2D 
connection of two manifolds. In this case, in the above shooting approach, we arrive at 
two equations for two parameters given in fl3.10p . In the analytic case p = 3, 5, where 
the dependence on parameters is assumed to be analytic, we conclude that the set of VSS 
profiles is at most countable. Moreover, for fixed p < Pq and any a > 0, we expect a finite 
(but, possibly, arbitrarily large for some values) number of VSS profiles according to the 
above p- and a-bifurcations and extra possible a-saddle-node bifurcations that much less 
are known about. 

On oscillatory blow-up in the ODE. Here we briefiy discuss an important aspect 
of ODEs such as fll.4p . = 1 (or in radial geometry). Namely, these admit a strong 
(nonlinear) unstability of orbits V{y) that is associated not with the exponential growth 
in the linear setting, but with the blow-up due to the main two terms in the ODE in (13. ip . 
For m = 2 these are 

(3.11) = -\V\P-^V, 
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Oscillatory blow-up in log-scale 




Figure 15. Oscillatory blow-up for the ODE (f3TTT]) for p = | and p = 2. 



where we omit all linear terms that are negligible as \V\ — > +00. Obviously, fl3.11l) does 
not admit blow-up solutions of constant sign, so we need to describe the oscillatory blow- 
up. In Figure [T5l we present typical blow-up behaviour of solutions In of (13. lip for 
y > with Cauchy data 



ViO) 



V'{0) = V"{0) = V"'{0) = 0. 



Solutions blow-up at some finite yo > with definite oscillatory behaviour as y — > ?/o- 

In order to detect this oscillatory behaviour, by shifting yo — y ^ —y, we have that 
blow-up happens at y = 0, i.e., as ?/ — > 0~, and introduce the oscillatory component ip{s) 
as follows: 

(3.12) yiy) = {—y)^fis), where s = \n{—y) and /i = — ^ < 0. 

Then ip{s) solves the following autonomous ODE in M: 



(3.13) 



+ 2(2^ - 3)^'" + (6;u2 - 18/i + ll)y?" + 2{2fi^ - V 
+ 11/i - 3)ip' + fiifi - l)(/i - 2)(/i - 3)ip 



Obviously, this ODE also admits blow-up singularities as ( 13. lip , but we are interested in 
bounded or slower growing orbits that can be extended up to the blow-up point s = — 00. 

In addition, using a standard min-max scenario of orbits behaviour on parameters, 
we expect a bounded (hopefully, periodic) or a spiral-type solution v^*(s) in M to exist, 
which, according to (13.121) . will describe a generic structure of blow-up singularities. The 
blow-up instabilities make it very difficult to see necessary orbits of equation (13.131) even 
numerically. A logarithmic-type distribution of zeros close to blow-up points as in (13.121) 
is seen in Figure [151 
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What is crucial in this case, is that the whole unstable bundle is two-dimensional 
incorporating the parameter yo as the position of the blow-up point (i.e., y ^ yo — y 
in (13.121) ) and the parameter sq of translation, ip^:{s) ^— ip^{s + sq). Therefore, our VSS 
profiles are those that have nothing to do with this singular blow-up bundle, so we again 
observe two equations for two parameters in (13.101) . 

4. Passing to the limit a ^ oo 

Without loss of generality, we consider the ID case in (11.41) and bear in mind the first 
symmetric similarity profile Vo{y) for the ODE 

(4.1) (-l)™+i1/(2") + ^yV^' + /5r- l^lf-V = 0, where f] = ^ ^ +oo 

as a — > +CXD. For /3 ^ 1, the first natural step is to use the scaling (13. 7p and to write the 
resulting equation in the perturbed form, 

(4.2) B^{V) = (-l)-+V(2-) + V- \V\P~'V = -| 2^ yV. 

Hence, the right-hand side becomes negligible as /3 ^ oo on smooth solutions. 
Formally passing to the limit (3 oo in (14.21) . we obtain the limit equation 

(4.3) Doo(Vr) = (-l)™+il^(2H + W - \W\P-^W = for y > 0, W{+oo) = 0. 
Let us show that a nontrivial solution of (14.31) does not exist. For instance, for m = 2, 

_-^y(4) + w- \W\P-^W = 0, so that, for y > 1, 

(4.4) 

= W+... =^ 3 ID bundle, W{y) ~ Ce-^ + ... . 
This ID is not enough to satisfy two symmetry conditions (13.41) . Similarly, for m = 3, 

W(^') + W -\W\P~^W = 0, andfor?/>l. 



(4.5) 



1^(6) = -w + ... =^ 3 2D bundle with char, values A± = -| ± i ^. 



The 2D bundle is not sufficient to satisfy three symmetry conditions at the origin, 

W'{0) = W"'{0) = W^^\0) = 0. 

The same lack of dimension happens for any m > 1. In other words, nonexistence for 
(14. 3 p is a typical property. 

Such nonexistence is easily checked numerically. In Figure Uni we show for m = p = 2 
that as a increases the VSS profiles with a fixed length of the flat part (such solutions 
exist) become more and more oscillatory and ceases to exist for a = +00. 

On the other hand, profiles with huge fiat non-oscillatory parts get very wide for a ^ 1. 
For instance. Figure [17] shows for p = 1.5 that for a = 40, the last less oscillatory VSS 
profile V2k{y) has the fiat part of the length ~ 380, while the first profiles are strongly 
oscillatory (and do not admit passage to the limit a 00). Notice that according to 
bifurcation exponents (11.70 . for a = 40, we expect, at least, 

164 profiles V2k and 164 profiles K-^, 
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Solutions of limit equation for large a, m=p=2 
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Figure 16. Towards nonexistence for the ODE (j4.4p . m = p = 2: solutions of 
(j4.2p get very oscillatory for a ^ 1. 
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Figure 17. Some VSS profiles satisfying (|42]) for m = 2, p = 1.5, and a = 40. 

i.e., 328 VSS profiles overall (and more in view of possible saddle- node bifurcations be- 
forehand, at some as~n ^ (0,40)). 

Therefore, as a — +oo, one can that the VSS profiles get wider indefinitely in the sense 
that, in the rescaled variables, uniformly on bounded intervals. 



(4.6) 



Vo{y) 



as a 



-oo, 



and the divergence is sufficiently (exponentially) fast. 

As a simple illustration of this behaviour presented in Figure [T71 we give a nonexistence 
result concerning a step-like profile for the limit equation (14.41) . Assume that, unlike 
(14. 6p . there exists a nontrivial finite limit as a — oo, so (14.41) admits solution W{y) with 

21 



(4.9) 



exponential decay such that 

(4.7) W{y)^l for ye[0,L], with some L > 1, 

and that the internal layer at ?/ ~ L does not depend essentially on L (recall that (14.41) is 
an autonomous ODE), and, in a natural sense, 

(4.8) J^W^<^L and \W\p+' L. 

The following nonexistence result has a "conditional" nature, where we prohibit certain 
solutions of a specific spatial shape: 

Proposition 4.1. The ODE i \4-4^ does not admit a nontrivial solution with exponential 
decay at infinity satisfying ( |^. 7[ ) and ^4-^ - 

Proof. We use two identities obtained by multiplying f l4.4p in by W and yW (Po- 
hozaev's multiplier in elliptic theory, [28]) to get 

-lj{w"f-\jw^ + ^j\wr' = ^. 

Substituting j {W")'^ from the first into the second identity yields 
(4.10) jW' = ^,j\Wr\ where 1, = 

Since 70 < 1 for any p > 1, this contradicts hypotheses (14.71) . (14. 8 p under which both 
integrals in (I4.10p are close to L. □ 

5. On some VSSs results for the non-monotone absorption 
In this short section, we orient our nonexistence business to the following PDE: 

(5.1) = -(-A)'"M-t"|M|f, 

where, unlike (11. ip . we take the non-monotone absorption nonlinearity \u\^- 

5.1. On local existence close to bifurcation points. Local existence of similarity 
solutions for p < Po for equation (11.40 . where 

(5.2) \V\P-^V ^\V\P, 

can be established as above by bifurcation theory. Then, for several types of "positively 
dominant" VSS profiles V{y) as in Figures [Tf [9l [TOl [T2l and[T3](a). their structure changes 
slightly after the replacement (15.20 . 

For instance, in Figure [181 "we present even VSS profiles for the ODE 

(5.3) -V^^^ + \yV' + (3V -\V\P = in M (/5=i±f), 

in the case p = 1.5, a = 7, where most of the solutions are very close to those in 
Figure [T3](a). Excluding first sufficiently small profiles, all the others are clearly positively 
dominant (the last ones look "almost" positive; careful checking shows that all of them 
remain oscillatory). Small solutions (dotted lines as well as that of order ~ 2 ■ 10~^ given 
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VSS profiles V with -|V|P, m=2, p=1 .5, a=7 
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Figure 18. Twenty nine VSS solutions V2k{y) of ([5l3]) for m = 2, p = 1.5, a = 7. 

VSS profiles V witfi -|V|P, m=2, p=1 .5, a=7 
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Figure 19. Seventeen VSS solutions V(j^{y) of (j5.3p for m = 2, p = 1.5, and a = 7. 

by the solid line) become essentially different. Instead of (13. 9p . equation (15. 3p admits the 
single constant non-zero equilibrium 

(5.4) V+ = +1 (and Vo = 0), 

so that VSS profiles K-^ such as in Figure [TT] and [12] are impossible. Since (15. 4p means 
that for the nonlinearity — two equilibria —1 and just coincide (in fact, identified), 
we obtain another but similar family of profiles again denoted by {V^i} that are shown in 
Figure [T9l For large a = 40, step-like solutions of (15. 3p are practically indistinguishable 
from those in Figure [T71 so that this behaviour as a — +oo suits both monotone and 
non-monotone absorption terms. 
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On the contrary, nonexistence of VSSs for 

(5.5) P>Po = l + 

can be proved much easier than for the original PDE (11. ip . This is explained in detail in 
[17] for a = 0, so we present here a few nonexistence comments. 

5.2. On nonexistence of similarity profiles. The similarity solutions remain the same, 
(11. 3p . with the elliptic equation f ll.4p . where we use the replacement (15. 2p . Therefore, 
integrating this equation over with the condition of exponential decay, we obtain in 
particular that 

(5.6) J\V\P = for p = po =^ V = 0. 

Therefore, all the continuous p-branches do not cross the vertical line p = Po, so that in 
the parameter region {p > pq} there are no VSS profiles at all. 
Obviously, for the original equation (11.41) this analysis leads to 

(5.7) J\V\P-^V = for p = po, 

which is not that controversial. Of course, this prohibits clearly positively dominant VSS 
profiles such as those in Figure fW a). but not obvious for profiles having essential negative 
parts as in (b) therein. Nevertheless, one can see that, for all these profiles, (15.71) is valid, 
since these are also positively dominant. 

5.3. Nonexistence for the PDE. Using a slight modification of the statement and the 
proof in [T71 § 2] based on Pohozaev's nonlinear capacity approach (see details and related 
references therein), we formulate the following nonexistence result for (15.10 : 

Theorem 5.1. Let a > —1 and 

(5.8) P>Po = l + 

// a function u{x, t) G ^^^^(M^ x R_|_ \ {(0, 0)}) satisfies (15. ip and (I in the weak sense, 
then M = 0. 

In particular, this means that, for this nonlinear ity, any similarity problem such as (11.41) 
does not admit a nontrivial solution. Of course, this implies that no other, non-similarity, 
VSSs exist. For the monotone nonlinearity as in (II. 4p . the proof of nonexistence gets more 
involved, and can be done, with some technical changes, along the lines in [TTl § 3.2]. 

Acknowledgement. The author thanks A.E. Shishkov for drawing his attention to 
time dependent absorption phenomena. 
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